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The ubiquitous occurrence of glycerol dialkyl glycerol tetraethers (GDGTs) in soils and 
their ability to record temperature and environmental changes offer the prospect of 
independently reconstructing continental paleotemperature and paleoenvironment from 
the loess-paleosol sequences (LPS) from the Chinese Loess Plateau (CLP). In this study 
we present records of GDGT-derived proxies for the last 70 kyr from the Yuanbao LPS, 
western CLP Temperature record reconstructed from the cyclization and methylation index 
of branched tetraethers (MBT-CBT) displays that the onset of deglacial warming at ~20 kyr 
before present (BP) precedes the strengthening of summer monsoon at ~ 1 5 kyr BP 
which is in agreement in timing with previous MBT-CBT temperature records from the 
southeastern CLP The maximal deglacial warming of ~10°C is slightly higher than those 
in the southeastern CLP perhaps due to the higher latitude and farther inland of the study 
site. The Branched and Isoprenoid Tetraether (BIT) index shows higher values (0.87-0.96 
range, 0.93 average) in the glacial loess and lower values (0.76-0.91 range, 0.83 average) 
in the Holocene paleosols, with a steady decreasing trend since the early Holocene. The 
decreasing trend could suggest enhanced Thaumarchaeota relative to GDGT producing 
bacteria activity since the early Holocene, but other possibilities, such as preferential 
degradation of isoprenoid GDGTs or upward increase in living archaea relative to bacteria 
in the paleosol profile, cannot be fully excluded. Our results thus demonstrate the need of 
future study on microbial community structure in soil column and differential degradation 
of GDGT molecules. 
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INTRODUCTION 

The microbial membrane lipids glycerol dialkyl glycerol 
tetraethers (GDGTs) are widespread, but different in structure, 
in archaea and bacteria (Figure 1). Archaea have isoprenoidal 
hydrocarbon chains, containing 0-3 cyclopentyl moieties or 
a cyclohexyl moiety in addition to four cyclopentyl moieties 
(i.e., crenarchaeol), at the sn 2,3 positions between the glycerol 
backbones (Schouten et al., 2000; Sinninghe Damste et al., 2002). 
Whereas some specific bacteria have branched alkyl chains with 
4-6 methyl groups, containing up to two cyclopentyl moieties, 
at sn 1,2 positions between the glycerol backbones (Weijers et al, 
2006a). 

Both the archaea and bacteria produced GDGTs have been 
used to develop temperature proxies that are promising and 
valuable in paleoclimatic studies. For example, sea surface tem- 
peratures were found being recorded by the distribution of iso- 
prenoidal GDGTs (i-GDGTs) of marine Thaumarchaeota with 
a varying number of clyclopentane moieties, and thus a proxy 
called TEX86 (TetraEther Index of tetraethers having 86 carbons) 
has been developed (Schouten et al., 2002). TEXs6 has also been 
used to estimate lake water temperatures in large lakes that receive 
minimal soil runoff (Blaga et al., 2009; Powers et al., 2010). Soon 
after, a proxy for continental air temperatures and soil pH was 
also proposed based on the distribution of soil bacterial branched 



GDGTs (br-GDGTs) by the definition of two indices, the methy- 
lation of br-GDGTs (MBT) and the cyclisation of br-GDGTs 
(CBT), as follows (Weijers et al., 2007): 

MBT = [la + lb + Ic]/[Ia + lb + Ic + Ha + lib + Ik + Ilia 

+ Illb + Ilk] (1) 
CBT = -LOG ([lb + IIb]/[Ia + Ila]) (2) 

The numbers in the equations refer to the abundance of the 
molecules in Figure 1 . The soil pH and mean annual air temper- 
ature (MAT) can be then estimated by the following relationship: 

CBT = 3.33 - 0.38 * pH (n = 114; r 2 = 0.70) (3) 
MBT = 0.122 + 0.187 x CBT + 0.020 

x MAT (n = 1 14; r 2 = 0.77) (4) 

Recently, the MBT index was adjusted according to a larger soil 
dataset (n = 278) as follows (Peterse et al., 2012): 

MBT' = [la + lb + Ic]/[Ia + lb + Ic + Ila + lib + lie + Ilia] 

(5) 
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FIGURE 1 | Molecular structures of isoprenoid GDGTs and branched GDGTs. 



in which compounds Illb and IIIc are not present, as in 
Equation (1), due to their undetectable amounts in most soil 
samples. As a result, a new calibration, generally in better 
agreement with independent proxy data, was proposed (Peterse 
etal.,2012): 



MAT' = 0.81 -5.67 x CBT + 31.0 x MBT' (n = 176; 



: 0.59) 
(6) 



The relative concentration of br-GDGTs and i-GDGTs in sed- 
iments and soils, potentially reflecting the relative br-GDGT 
producing bacteria and archaea input, is also a promising envi- 
ronmental indicator. Hopmans et al. (2004) introduced the 
Branched and Isoprenoid Tetraether (BIT) index for tracing soil 
organic matter (OM) input into a marine system as follows: 



BIT = [I + II + III] /[I + II + III + Cren] 



(7) 



where I, II, and III are the abundances of br-GDGTs without 
cyclopentane moieties and Cren is the abundance of crenarchaeol 
in a given sample (structures in Figure 1). The BIT index there- 
fore ranges from 0 to 1. A value of 0 indicates the presence of 
only crenarchaeol with no br-GDGT inputs, and a value of 1 indi- 
cates the presence of only br-GDGTs, which was thought initially 
to represent marine OM and soil OM, respectively (Hopmans 
et al., 2004). However, subsequent works showed that crenar- 
chaeol can be produced in small to moderate amounts in soils 
(Leininger et al, 2006; Weijers et al, 2006a,b; Walsh et al, 2008) 
and that br-GDGTs may also be produced in situ in aquatic envi- 
ronments (Peterse et al, 2009; Sinninghe Damste et al., 2009; 
Tierney and Russell, 2009; Tierney et al., 2011; Zhu et al, 2011; 
Wang et al, 2012; Zhang et al., 2012). As an extreme example, in 
highly alkaline soils (pH > 7.5) under severe drought conditions, 
archaeal i-GDGTs with crenarchaeol as the main component were 



found predominant over bacterial br-GDGTs, as indicated by the 
increase in the abundance ratio of i-GDGTs to br-GDGTs, i.e., 
Ri/b, to values as high as 6 at higher pH values (Xie et al., 2012). 
So the BIT Index might be difficult to accurately trace soil OM 
in aquatic environments under drought conditions. Nevertheless, 
the BIT Index, as well as R;/b value, may potentially be used to 
trace variations in microbial structure in terms of archaea ver- 
sus br-GDGTs producing bacteria, as results of soil environmental 
changes. 

The ubiquitous occurrence of GDGTs in soils offers the 
prospect of independently reconstructing continental pale- 
otemperature from terrestrial deposits, e.g., the loess-paleosol 
sequences (LPS) from the Chinese Loess Plateau (CLP) 
(Figure 2). The LPS resulted from changes in monsoon inten- 
sity; loess is mainly deposited during cool and dry periods when 
East Asian winter monsoon (EAWM) intensifies, whereas soil 
development takes place during the warmer and wetter periods 
when East Asian summer monsoon (EASM) builds up (An, 2000). 
Although various climate proxies, e.g., pedogenic magnetic sus- 
ceptibility (MS), phytoliths, grain size distributions and S ls O 
of rhizoconcretions and land snail shells, have been applied to 
unravel the history of the East Asian monsoon for the past sev- 
eral million years (e.g., Maher et al., 1994; Xiao et al., 1995; Li 
et al., 2007; Lu et al., 2007), little is known about the temper- 
ature history in the CLP region until the successful application 
of br-GDGTs for temperature reconstruction on the southeast 
CLP recently (Peterse et al., 2011; Gao et al., 2012). However, 
these works are clearly insufficient to reveal the whole picture 
of the temporal and spatial paleotemperature distributions with 
regard to the vast area of the CLP. In addition, failure examples 
of applying MBT-CBT proxy to three well-studied LPS in Serbia, 
northeast Siberia, and on Mt. Kilimanjaro prompt us to be cau- 
tious when interpreting these indices in LPS studies (Zech et al., 
2012). 
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FIGURE 2 | Chinese Loess Plateau (CLP) and the site of Yuanbao (YB) 
in this study. Sites LT (Lantian) and MS (Mangshan) have been studied for 
MBT-CBT temperature records by others and are mentioned in this paper. 



In this study, we present GDGTs results for the past 70 kyr from 
the Yuanbao LPS, western CLP. In addition to MBT-CBT tem- 
perature reconstruction and its comparison to previous results, 
records of BIT and R;/b indices are also reported and discussed, 
which were not included in previous works. 

MATERIALS AND METHODS 
SITE AND SAMPLING 

The site of Yuanbao sequence (103°09'E, 35°38'N) is located near 
the Linxia city in the western edge of the CLP prevailing with a 
typical temperate semi-arid climate (Figure 2) (Chen et al., 1997). 
The MAT at the Linxia city is ~7.0°C, with winter temperature 
below freezing (DJF) (Figure 3). The precipitation is concen- 
trated from late spring to early autumn (May to September), 
accounting for 80% of the annual mean of 500 mm (Figure 3); 
whereas the annual evaporation is 1300 mm. The study site lies on 
the fourth terrace of the Xiahe River with an elevation of 2177 m, 
which is about 500 m above the Linxia city. According to the lapse 
rate of 0.6°C/100 m, the MAT at the site is estimated to be ~4°C. 
The loess on the river terrace is up to 150 m thick. 



In the present study, we investigated the upper 27 m of the 
Yuanbao sequence from the bottom of loess Li to the paleosol 
So with the top sample 0.5 m below the surface. The sequence 
was spliced by two exposed sections, CH02/02 and CH02/03, with 
2-3 km apart (Chen et al., 1997, 1999; Lai and Wintle, 2006; Lai 
et al., 2007). The two sections have been dated by using high reso- 
lution optically stimulated luminescence (OSL), showing that the 
4.0-m long section CH02/02 have been deposited since ~17.3 kyr 
before present (BP) and the interval from 17 m up to 0.5 m in 
the section CH02/03 was deposited from 51.2 to 17.5 kyr BP 
(Lai and Wintle, 2006; Lai et al., 2007). The exposed outer lay- 
ers of the sections were removed during the course of sampling 
to ensure fresh samples for this study. A total of 40 samples, 15 
in the section CH02/02 and 25 in the section CH02/03 were col- 
lected. The uppermost 0.5 m of in the CH02/02 and 1.1 m of the 
CH02/03 section are disturbed by ploughing, and thus were not 
sampled. 

MAGNETIC SUSCEPTIBILITY, BULK ORGANIC CARBON AND SOIL pH 
ANALYSIS 

Samples were oven dried at 50° C, lightly ground and aliquots 
of ~10g were analyzed using a Bartington MS2 for MS. Total 
organic carbon (C or g) was analyzed in duplicate using an elemen- 
tal analyzer (Vario Pyro Cube) after removal of carbonate with 
diluted HC1. The average standard deviations of C org is ±0.2%. 
The pH of the paleosol and loess samples was measured in a 1:2.5 
soikwater (w/v) mixture using a digital pH meter. 

GDGT ANALYSIS 

All air-dried samples, with each ca. 50 g in weight, were homog- 
enized with a mortar and pestle, and ultrasonically extracted 
(9 x) with MeOH (3 x), dichloromethane (DCM)/MeOH (1:1, 
v/v; 2 x ) and DCM (2 x ) and all extracts were combined after 
centrifugation. Known amounts of an internal C46 GDGT stan- 
dard were added according to Huguet et al. (2006). The total 
extracts were dried by way of rotary evaporation under vacuum 
and then were purified and separated into an apolar fraction 
and a polar fraction over an activated silica gel column by elu- 
tion with n-hexane and DCM/MeOH (1:1, v/v), respectively, the 
GDGTs being in the polar fraction. The solvent was removed 
under N2 and the residue dissolved via sonication (5 min) in hex- 
ane/propanol (99:1, v/v) and filtered through a 0.45 [im, 4 mm 
diameter PTFE filter. 

Analysis of GDGTs was performed using an Agilent 1200 
HPLC/6410 TripleQuad MS instrument equipped with an 
auto-injector and Chemstation chromatography manager soft- 
ware. Procedures described by Schouten et al. (2007) were 
applied. Separation was achieved with a Prevail Cyano col- 
umn (2.1 x 150 mm, 3 (Jim diameter; Grace, USA), maintained 
at 30°C. Injection volume varied from 10 to 20 uT GDGTs 
were eluted isocratically with 99% hexane and 1% propanol 
for 5 min, followed by a linear gradient to 1.8% propanol in 
45 min. Flow rate was 0.2 ml/min. Detection was achieved using 
atmospheric pressure chemical ionization mass spectrometry 
(APCI-MS) via selected ion monotoring (SIM) of [M + H] + 
ions (in MSI) and GDGTs were quantified by integration 
of the peak areas. Absolute quantification was achieved by 
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FIGURE 3 | Monthly distributions of temperature and precipitation 




between 1971 and 2000 at the Linxia city. Data are from Public 




Meteorological Service Center, China Meteorological Administration. Note 


that the Yuanbao section is about 500 m above the Linxia city. 
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calculating the area of the corresponding peaks in the chro- 
matograms, comparing them with the peak area of the internal 
standard. 



RESULTS AND DISCUSSION 
AGE MODEL, MAGNETIC SUSCEPTIBILITY AND C 



org 



CONTENT 



Age determination was according to the OSL dating results on the 
sections CH02/02 and CH02/03 (Lai and Wintle, 2006; Lai et al, 
2007). For the section CH02/02, our MS curve was correlated 
with that from Lai and Wintle (2006) before age determination 
instead of applying the given depth-age relationship, which may 
reduce the errors from depth measurement difference between 
the two studies. For the section CH02/03, because no MS data 
were presented along with age determination in Lai et al. (2007), 
the depth-age relationship (Lai et al., 2007) was applied to our 
samples. 

Based on the age model, MS and C org (wt%) content in the 
Yuanbao sequence changed in parallel, with lower values dur- 
ing the glacial time, the deglacial onset of increases at ~15kyr 
BP, and higher values during the Holocene (Figures 4B,C). The 
two records are similar to the changes of speleothem 8 ls O 
(Figure 4A), as is expected because all of these proxies are asso- 
ciated with the development of the EASM precipitation (Zhou 



et al, 1990; Wang et al, 2001; Dykoski et al, 2005; Peterse et al, 
2011). However, there is an obvious asynchrony between the 
peaks of MS (at 2.2 kyr BP) and the C org content (at 6.9 kyr BP). 
This mismatch in time might be due to that MS is a complex 
climate proxy and its value is not always consistent with the 
soil development, which has been pointed out by Guo et al. 
(1998). The low sampling resolution might be also responsible 
for the mismatch. But overall, records of MS and C org suggest 
that EASM was strongest during the mid Holocene, consistent 
with numerous results on the regional paleoclimate [reviewed by 
(Zhao et al., 2009)]. A fine comparison of the timing of EASM 
development with other records, e.g. the Mangshan record, is dif- 
ficult at present, because of the low sampling resolution in this 
study and the different dating methods for age models. However, 
the onset of EASM enhancement seems synchronous across CLP, 
as both MS and C or g began to rise at ~ 15 kyr BP in both Yuanbao 
and Mangshan sequences. 

GDGTs CONCENTRATION AND DISTRIBUTION 

A remarkable feature of GDGT concentration profiles in the 
Yuanbao LPS is the large differences in their concentration and 
distribution between the last glacial loess and the Holocence 
paleosol (Supplementary data sheet). 




FIGURE 4 | Variations down the loess-paleosol sequence at Yuanbao 
section of (B) magnetic susceptibility, (C) C org content, (D) 
concentrations of i-GDGTs and br-GDGTs in dry samples, and (E) 
Com -content normalized concentrations of i-GDGTs and br-GDGTs. 



Speleothem 8 ls O (A) indicating East Asian summer monsoon (Wang et al., 
2001 ; Dykoski et al., 2005) is applied for comparison. The gray area denotes 
the Holocene period since ~12 kyr BP, and the area before that denotes the 
last glacial time. 
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Br-GDGTs are present in all paleosols, but they were detected 
in only 19 out of 31 loess samples. Their total concentration 
in the paleosols varies between 2.9 and 19.6 ng g -1 , whereas it 
is much lower in loess samples with values between 0.04 and 
0.9 ng g _1 (Figure 4D). However, if normalized to TOC content, 
the difference between loess and paleosol becomes reduced, i.e., 
0.5-3.5 (xg g _1 TOC in paleosols and 0.05-2.0 |xg g _1 TOC in 
loess (Figure 4E). In addition to the concentration difference, 
loess samples do not contain detectable amounts of all nine br- 
GDGTs that are used for the CBT and MBT calculations, as 
paleosols do. In particular, br-GDGTs Ic, lie, and Ilk in loess 
samples are either completely undetectable (i.e. br-GDGT Ilk) 
or frequently below detection limit and, if present, less than 1% 
of total br-GDGTs on average (i.e., br-GDGTs IC and IIC). The 
most dominant br-GDGTs in the LPS is Ha, accounting for 33.2% 
(23.4-45.2% range) in total br-GDGTs. Br-GDGT la, Ilia, lb and 
lib are present in all samples, but with smaller percentages (3.6- 
25.1% range). Br-GDGT Illb is less than 8.4% and in several loess 
samples it is below detection limit. 

TGDGTs also occur, although in less amounts, in the 
sequence. Their total concentrations are higher in paleosols (0.6- 
4.0ng g _1 dry bulk sample, or 0.09-0.76 |xgg _1 TOC) than in 
loess (0.001-0.06 ng g" 1 dry bulk sample, or 0.002-0.12 u,g g" 1 
TOC) (Figures 4D,E). All i-GDGTs in Figure 1 are present in 
all the paleosols with crenarchaeol the dominant component 
(40.5-65.7% range); however, crenarchaeol becomes the only i- 
GDGT compound in the loess except two samples. This suggests 
that archaeal community structure in the glacial loess might be 
quite different from that in the Holocene paleosols, with the 
crenarchaeol-producing Thaumarchaeota predominant during 
the last glaciation. The cause for the Thaumarchaeota predom- 
inance in archaeal community might be related to the slightly 
more alkaline conditions during the glacial time, as indicated 
by the higher pH values measured in loess samples (Figure 5F), 
under which elevated ammonia oxidation rate may promote the 
growth rate of soil ammonia-oxidizing Thaumarchaeota (Xie 
et al.,2012). 

The occurrence of higher br-GDGT concentration in pale- 
osol than in loess has been also observed in the Mangshan 
sequence (Peterse et al., 2011), although not observed in the 
Lantian sequence (Gao et al., 2012). Higher concentrations of 
GDGTs, including both i-GDGTs and br-GDGTs, in paleosol sug- 
gest that microbial activity is an important agent in the pedogenic 
processes. Likewise, the higher TOC normalized concentrations 
of GDGTs in paleosols could indicate greater contributions of 
microbial OM to the total OM. However, although GDGTs were 
not detected in some loess samples in this study, they were present 
in the loess layers at the Mangshan and Lantian sequences (Peterse 
et al., 2011; Gao et al., 2012). Two factors during the last glacial 
time might be responsible for the discrepancy: (1) Temperature 
was lower, and perhaps aridity was higher, at Yuanbao than at 
Mangshan and Lantian, which was more unfavorable for micro- 
bial and biomass growth. (2) Dust accumulation rate was several 
time higher in the NW CLP than in the SE CLP (Lai et al, 2007). 
So GDGTs concentrations could have been more diluted in the 
Yuanbao loess layers, and thus are hard to be detected unless loess 
sample amounts are added. 



TEMPERATURE RECORD 

Although atmospheric dust from the loess plateau has not been 
analyzed, samples in loess source regions and a dust sample from 
near the west coast of Central Africa have been found containing 
negligible amounts of GDGTs (Hopmans et al., 2004; Gao et al., 
2012). Thus, GDGTs in the LPS are likely locally produced, and 
the MBT-CBT record reflects local conditions (Peterse et al., 201 1; 
Gao et al.,2012). 

Due to lack of surface soil sample, temperature estimate of 
~8.3°C from Equation (4) or 7.6° C from Equation (6) for the 
youngest sample (i.e., ~0.9kyr BP) is compared to present. It 
is higher than present MAT of ~4°C, but lower than summer 
(June- August) mean of ~14°C at the Yuanbao site. In the Lantian 
sequence, temperature estimate of ~ 18.4° C from Equation (4) for 
the youngest sample at 2.9 kyr BP also fell in between the local 
MAT of 13.1°C and summer mean of 26°C (Gao et al, 2012). In 
the Mangshan sequence, estimate of ~23°C from Equation (4) 
for samples for the last 2 kyr was in agreement with local sum- 
mer temperature (Peterse et al., 2011). Reconstructed MATs bias 
to local warm season temperatures have been also observed in 
some investigations on sediments from lakes, coastal area and 
peat bog (Rueda et al, 2009; Huguet et al, 2010; Sun et al, 201 1). 
Thus it is likely that the GDGT-producing bacteria might thrive in 
the warm season of a year. However, in a 1 -year investigation on 
soils sampled from eight different plots in the mid-latitude USA, 
The Netherlands and the UK, the MBT-CBT temperature proxy 
did not show any seasonal patterns due to the fact that branched 
GDGT core lipids represent a standing stock that has accumulated 
over the course of years (ca. 20 years, Weijers et al, 201 1). In that 
investigation, the authors also found that the MBT-CBT temper- 
ature, though calibrated with annual MAT, did not always exactly 
reflect the MAT, which resulted from the offset between soil and 
air temperature (Weijers et al., 2011). Therefore, the higher than 
global mean offset between soil and air temperature in the CLP 
region might be the cause of higher MBT-CBT temperature than 
the local MAT observed here and in previous works (Peterse et al., 
2011, 2012; Gao et al., 2012). Nevertheless, the MBT-CBT proxy 
is still thought to be able to provide reasonable estimates of past 
changes of MAT, although the reconstructed absolute tempera- 
ture is associated with a slightly larger error (ca. 5°C) (Weijers 
etal., 2011). 

The MBT-CBT derived temperature calculated by application 
of Equation (4) for the Yuanbao sequence varies between 1.1°C 
and 20.8°C. Temperature was <8.1°C in the last glacial time 
between 70 and 20 kyr BP, and then a deglacial warming occurred 
from 3.7°C at 19.7 kyr BP to 20.8°C at 14.3 kyr BP, the highest 
value in the time frame (Figure 5C). The Holocene is featured 
by a continuous cooling trend towards a temperature value of 
8.3°C for the top sample. Since br-GDGT Ilk and IIIB are either 
undetectable or a minor component in our samples, Equation (5) 
was also used to recalculate the MBT-CBT derived temperature. 
The air temperatures that are derived with the new calibration are 
lower than the original record for most samples, but higher for 
several loess samples in the last glacial time (Figure 5C). Thus, the 
new record suggests a smaller increase in air temperature over the 
deglaciation with temperature estimates increasing from 5.1°C to 
a maximum of 15.8°C. 
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FIGURE 5 | GDGT derived proxy records down the sequence at Yuanbao 
section and their comparison to speleothem 8 ls O. (A) Speleothem 8 18 0 
from Dongge and Hulu cave (Wang et al., 2001 ; Dykoski et al., 2005); (B) 
magnetic susceptibility; (C) MBT-CBT derived temperature. Black and gray 



curves are from Yuanbao sequence using Equations (6) and (4), respectively; 
red and blue curves are from Manshan (Peterse et al., 2012) and Lantian 
sequences (Gao et al., 2012), respectively; (D) BIT; (E) R i/b ; (F) Soil pH. 
OM; onset of East Asian summer monsoom; OW, onset of warming. 



The MBT-CBT derived temperature records from other LPSs 
also suggested a large increase in air temperature over the 
deglaciation, e.g., from 13.5°C at 20.5 kyr BP to a maximum 
of 27.2°C at 12.4 kyr BP in the Mangshan sequence and from 
14.6°C at 22.5 kyr BP to a maximum of 23.3°C at 12.2 kyr BP in 
the Lantian sequence (Peterse et al., 2011; Gao et al., 2012). By 
using the new calibration (6), Peterse et al. (2012) gave a smaller 
deglacial warming of 6 ~ 7°C for the Mangshan sequence (Red 
curve in Figure 5C). Our temperature record shows a deglacial 
increase of ~18°C or ~10°C by application of Equation (4) or 
(6), respectively. The deglacial increase of ~10°C from Equation 
(6) seems more reasonable, as 4.5-9°C cooler air temperatures 
during the LGM and 1-2° C warmer during the Holocene cli- 
matic optimum than at present in central Asia have been esti- 
mated based on pollen and phytolith studies (Sun et al., 1997; 
Lu et al, 2007). The lower absolute temperature estimates and 
larger deglacial warming at Yuanbao than those at Mangshan and 
Lantian may be due to the higher latitude and farther inland of 
the Yuanbao site. 

Similar to the findings of early warming on the SE CLP rela- 
tive to EASM precipitation increase by Peterse et al. (2011) and 
Gao et al. (2012), our record also suggests early warming in the 



western CLP, as indicated by the lead of onset of warming at 
~20kyr BP to the onset of MS and C org increase at ~15ky BP 
from the same sequence (Figures 5 A-C). The high temperature 
and still low EASM precipitation during the early Holocene is in 
agreement with the finding that wind-blown sands were mobile 
at many sites along the desert margin in northern China due to 
hot and dry conditions (Mason et al, 2009). Therefore, although 
the EASM precipitation history on its temporal and spatial distri- 
butions are still in debate (An et al, 2000; Zhao et al, 2009), the 
early deglacial warming was likely regionally synchronous across 
the CLP. 



Br-GDGTs vs. i-GDGTs 

The BIT index exhibits higher values (0.87-0.96 range, 0.93 aver- 
age) in the glacial loess and lower values (0.76-0.91 range, 0.83 
average) in the Holocene paleosols (Figure 5D). Similar occur- 
rence, but in opposite direction, is also shown in the record of 
R;/b that indicates the relative abundance of total i-GDGTs to total 
br-GDGTs (Figure 5E). As the production of both archaea andbr- 
GDGT producing bacteria could have increased in the pedogenic 
process during the Holocene, the records of BIT and R;/b could 
suggest that the increase of archaea might have exceeded that of 
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the br-GDGT producing bacteria. The increase of R;/b has been 
proposed to indicate alkaline soils induced by enhanced aridity 
(Xie et al., 2012), perhaps due to that soil Thaumarchaeota may 
produce higher amounts of crenarchaeol under alkaline condi- 
tions (Weijers et al, 2006b; Pearson et al, 2008). Indeed, CBT- 
derived pH record shows averagely higher values (~8.5 ±0.1) 
during the Holocene than during the last glciation (~8.1 ± 0.4) 
(Figure 5F). However, the pH record is not consistent with what 
was obtained from direct pH measurement on samples, showing 
slightly higher pH during the last glaciation. This disagreement is 
not clear at present, but might reflect the difficulty of the CBT 
index to precisely estimate alkaline soil pH (Xie et al., 2012). 
In addition, the increasing trend in R;/b record since the early 
Holocene is not in accordance with the unchanged CBT-derived 
pH record, or even in opposite direction to the directly measured 
pH record. Therefore, the Ri/b record in this study seems unlikely 
to be caused by soil pH changes. 

We notice that the steadily increasing (or increasing) trend 
of R;/b (or BIT) since the early Holocene matches well with the 
MBT-CBT temperature record (Figure 5C). However, it seems 
unlikely to establish a causal relationship between them, because 
other soil physicochemical conditions in addition to temperature, 
such as moisture, pH, redox conditions and nutrient levels, might 
be also important to influence the microbial community. Soil 
physicochemical conditions should have been associated more 
closely with precipitation and OM accumulation, which can be 
traced by MS and TOC records in this study. However, the BIT 
and Ri/b records for the Holocene do not match well with those 
of MS and C 0 i- g . 

Alternatively, the decreasing (or increasing) trend of BIT (or 
R;/b) since the early Holocene may be caused by the differential 
degradation of i-GDGTs vs. br-GDGTs. If so, the trend would sug- 
gest that i-GDGTs were preferentially degraded and/or br-GDGTs 
were selectively preserved in the sediment sequence. This infer- 
ence is in agreement with observations from marine sediment 
analysis (Huguet et al., 2008, 2009) and modeling experiments 
(Ding et al., 2013). For example, Huguet et al. (2008) indi- 
cated that degradation rates of crenarchaeol are 2-fold higher 
than those of soil-derived br-GDGTs. Besides, there are also 
other potential unrealized occurrences associated with the sec- 
ular BIT and Ri/b changes. For example, we do not know the 



depth distribution of living archaea and GDGT-producing bac- 
teria in soils, i.e., the growth depth effect as proposed by Zech 
et al. (2012). Therefore, upward increase in living archaea relative 
to bacteria in the paleosol profile cannot be fully excluded. This 
issue may be resolved by analyzing intact polar GDGTs, markers 
for living cells, in the profile in future. In any case, it reminds us 
to be cautious with BIT and R;/b proxies in paleo-environmental 
reconstructions. 

CONCLUSIONS 

The Yuanbao loess-paleosol sequence from the west CLP shows 
clear biogeochemical and climate contrasts between the last 
glaciation and the Holocene. Higher C org content and GDGT 
concentrations in the Holocene paleosols consistently indicate an 
intensified pedogenic process as results of strengthened warm and 
wet EASM. MBT-CBT derived temperature record displays an 
early deglacial warming relative to the onset of EASM, which is 
in agreement with previous findings from the southeastern CLP 
and suggests that the early deglacial warming was likely region- 
ally synchronous across the CLP. In the BIT (or Ri/b) variations, 
a steady decrease (or increase) since the early Holocene is promi- 
nent. Changes in microbial community, i.e., archaea production 
exceeding that of the br-GDGT producing bacteria, could be 
responsible for this. However, other possibilities, such as prefer- 
ential degradation of isoprenoid GDGTs or upward increase in 
living archaea relative to bacteria in the paleosol profile, cannot 
be fully excluded. Thus, future studies, with foci on microbial 
community structure in soil column and differential degradation 
of GDGT molecules, are needed for accurate paleoclimate and 
paleo-biogeochemical reconstructions from the CLP. 

ACKNOWLEDGMENTS 

We thank Drs. liawu Zhang and Suping Wang for collaboration in 
collecting loess and soil samples. This work is partly supported by 
NSF of China (Grant No. 41276072). Three reviewers are thanked 
for their valuable comments. 

SUPPLEMENTARY MATERIAL 

The Supplementary Material for this article can be found online 
at: http://www.frontiersin.org/Terrestrial_Microbiology/10.3389/ 
fmicb.20 13.001 99/abstract 



REFERENCES 

An, Z. (2000). The history and variabil- 
ity of the East Asian paleomonsoon 
climate. Q. Sci. Rev. 19, 171-187. 
doi: 10.1029/2010JG001365 

An, Z., Porter, S. C, Kutzbach, J. E., 
Wu, X., Wang, S., Liu, X., et al. 
(2000). Asynchronous Holocene 
optimum of the East Asian mon- 
soon. Q. Sci. Rev. 19, 743-762. doi: 
10.1016/S0277-3791(99)00031-1 

Blaga, C. I., Reichart, G.-J., Heiri, O., 
and Sinninghe Damste, J. S. (2009). 
Tetraether membrane lipid distri- 
butions in water-column particulate 
matter and sediments: a study of 47 
European lakes along a north-south 



transect. /. Paleolimnol. 41, 
523-540. doi: 10.1007/sl0933-008- 
9242-2 

Chen, F. H., Bloemendal, J., Feng, Z. 
D., Wang, J. M., Parker, E., and Guo, 
Z. T. (1999). East Asian monsoon 
variations during Oxygen Isotope 
Stage 5: evidence from the north- 
western margin of the Chinese loess 
plateau. Q. Sci. Rev. 18, 1127-1135. 
doi: 10.1016/S0277-3791{98) 
00047-X 

Chen, F. H., Bloemendal, J., Wang, 
J. M., Li, J. J., and Oldfield, F. 
(1997). High-resolution multi- 
proxy climate records from 
Chinese loess, evidence for rapid 



climatic changes over the last 75 
kyr. Palaeogeogr. Palaeoclimatol. 
Palaeoecol. 130, 323-335. doi: 
10.1016/S0031-0182{96)00149-6 

Ding, W. H., Yang, H., He, G. Q., and 
Xie, S. C. (2013). Effect of oxidative 
degradation by hydrogen peroxide 
on tetraethers based organic prox- 
ies. Q. Sci 33, 39-47. 

Dykoski, C. A., Edwards, R. L., Cheng, 
H., Yuan, D., Cai, Y, Zhang, M., 
et al. (2005). A high-resolution, 
absolute- dated Holocene and 
deglacial Asian monsoon record 
from Dongge Cave, China. Earth 
Planet. Sci. Lett. 233, 71-86. doi: 
10.1016/j.epsl.2005.01.036 



Gao, L., Nie, J., Clemens, S., Liu, W., 
Sun, J., Zech, R., et al. (2012). The 
importance of solar insolation on 
the temperature variations for the 
past 110 kyr on the Chinese Loess 
Plateau. Palaeogeogr. Palaeoclimatol. 
Palaeoecol. 317-318, 128-133. 
doi: 10.1016/j.palaeo.2011. 
12.021 

Guo, Z., Liu, T., Fedoroff, N., Wei, 
L., Ding, Z., Wu, N., et al. (1998). 
Climate extremes in loess of China 
coupled with the strength of deep- 
water formation in the North 
Atlantic. Glob. Planet. Change 18, 
113-128. doi: 10.1016/S0921-8181 
(98)00010-1 



www.frontiersin.org 



July 2013 | Volume 4 | Article 199 | 7 



Jia et al. 



Tetraether records in loess-paleosol sequence 



Hopmans, E. C.,Weijers, J. W. H., 
Schefuss, E., Herfort, L., Sinninghe 
Damste, J. S., and Schouten, S. 
(2004). A novel proxy for ter- 
restrial organic matter in sedi- 
ments based on branched and 
isoprenoid tetraether lipids. Earth 
Planet. Sci. Lett. 224, 107-116. doi: 
10.1016/j.epsl.2004.05.012 

Huguet, C, de Lange, G. J., Gustafsson, 
O., Middelburg, J. J., Sinninghe 
Damste, J. S., and Schouten, 
S. (2008). Selective preserva- 
tion of soil organic matter 
in oxidized marine sediments 
(Madeira Abyssal Plain). Geochim. 
Cosmochim. Acta 72, 6061-6068. 
doi: 10.1016/j.gca.2008.09.021 

Huguet, A., Fosse, C, Laggoun- 
Defarge, E, Toussaint, M.-L., and 
Derenne, S. (2010). Occurrence 
and distribution of glycerol dialkyl 
glycerol tetraethers in a French peat 
bog. Org. Geochem. 41, 559-572. 
doi: 10.1016/j.orggeochem.2010. 
02.015 

Huguet, C, Hopmans, E. C, Febo- 
Ayala, W., Thompson, D. H., 
Sinninghe Damste, J. S., and 
Schouten, S. (2006). An improved 
method to determine the absolute 
abundance of glycerol dibiphy- 
tanyl glycerol tetraether lipids. 
Org. Geochem. 37, 1036-1041. doi: 
10.1016/j.orggeochem.2006.05.008 

Huguet, C, Kim, J. H., de Lange, 
G. J., Sinninghe Damste, J. 
S., and Schouten, S. (2009). 
Effects of long term oxic degra- 
dation on the U-37(K '), TEX86 
and BIT organic proxies. Org. 
Geochem. 40, 1 188-1 194. doi: 
10.1016/j.orggeochem.2009.09.003 

Lai, Z. P., and Wintle, A. G. (2006). 
Locating the boundary between the 
Pleistocene and the Holocene in 
Chinese loess using luminescence. 
Holocene 16, 893-899. doi: 10.1191/ 
0959683606hol980rr 

Lai, Z. P., Wintle, A. G., and Thomas, 
D. S. G. (2007). Rates of dust 
deposition between 50 ka and 
20 ka revealed by OSL dating at 
Yuanbao on the Chinese loess 
plateau. Palaeogeogr. Palaeoclimatol. 
Palaeoecol. 248, 431-439. doi: 
10.1016/j.palaeo.2006.12.013 

Leininger, S., Urich, T, Schloter, M., 
Schwark, L., Qi, J., Nicol, G. W., 
et al. (2006). Archaea predom- 
inate among ammonia oxidizing 
prokaryotes in soils. Nature 442, 
806-809. doi: 10.1038/nature04983 

Li, G., Sheng, X., Chen, J., Yang, J., 
and Chen, Y. (2007). Oxygen- 
isotope record of paleorainwater in 
authigenic carbonates of Chinese 
loess-paleosol sequences and 
its paleoclimatic significance. 



Palaeogeogr. Palaeoclimatol. 
Palaeoecol. 245, 551-559. doi: 
10.1016/j.palaeo.2006.10.003 

Lu, H.-Y, Wu, N.-Q., Liu, K.-B., 
Jiang, H., and Liu, T.-S. (2007). 
Phytoliths as quantitative indi- 
cators for the reconstruction of 
past environmental conditions 
in China II: palaeoevironmental 
reconstruction in the Loess Plateau. 
Q. Set. Rev. 26, 759-772. doi: 
10.1016/j.quascirev.2006.10.006 

Maher, B. A., Thompson, R., and Zhou, 
L. P. (1994). Spatial and tempo- 
ral reconstructions of changes in 
the Asian palaeomonsoon: a new 
mineral magnetic approach. Earth 
Planet. Sci Lett. 125, 461-471. 
doi: 10.1016/0012-821X(94) 
90232-1 

Mason, J. A., Lu, H., Zhou, Y, Miao, 
X., Swinehart, J. B., Liu, Z., et al. 
(2009). Dune mobility and aridity 
at the desert margin of northern 
China at a time of peak monsoon 
strength. Geology 37, 947-950. doi: 
10.1130/G30240A.1 

Pearson, A., Pi, Y, Zhao, W., Li, 
W., Li, Y, Inskeep, W., et al. 

(2008) . Factors controlling the dis- 
tribution of archaeal tetraethers 
in terrestrial hot springs. Appl. 
Environ. Microbiol. 74, 3523-3532. 
doi: 10.1 128/ AEM.02450-07 

Peterse, E, Kim, J. H., Schouten, 
S., Kristensen, D. K., Koc, N., 
and Sinninghe Damste, J. S. 

(2009) . Constraints on the 
application of the MBT/CBT 
palaeothermometer at high latitude 
environments (Svalbard, Norway). 
Org. Geochem. 40, 692-699. doi: 
10.1016/j.orggeochem.2009.03.004 

Peterse, E, Prins, M. A., Beets, C. J., 
Troelstra, S. R., Zheng, H., Gu, Z., 
et al. (2011). Decoupled warming 
and monsoon precipitation in East 
Asia over the last deglaciation. Earth 
Planet. Sci. Lett. 301, 256-264. doi: 
10.1016/j.epsl.2010.1 1.010 

Peterse, E, van der Meer, J., Schouten, 
S., Weijers, J. W. H., Fierer, N., 
Jackson, R. B., et al. (2012). Revised 
calibration of the MBT-CBT 
paleotemperature proxy based on 
branched tetraether membrane 
lipids in surface soils. Geochim. 
Cosmochim. Acta 96, 215-229. doi: 
10.1016/j.gca.2012.08.011 

Powers, L., Werne, J. P., Vanderwoude, 
A. J., Sinninghe Damste, J. S., 
Hopmans, E. C, and Schouten, 
S. (2010). Applicability and 
calibration of the TEXg6 pale- 
othermometer in lakes. Org. 
Geochem. 4 1 , 404-4 13. doi: 
10. 10 16/j.orggeochem.2009. 1 1 .009 

Rueda, G., Rosell-Mele, A., Escala, M., 
Gyllencreutz, R., and Backman, J. 



(2009). Comparison of instrumen- 
tal and GDGT-based estimates 
of sea surface and air tem- 
peratures from the Skagerrak. 
Org. Geochem. 40, 287-291. doi: 
10.1016/j.orggeochem.2008. 10.012 

Schouten, S., Hopmans, E. C, Pancost, 
R. D., and Sinninghe Damste, 
J. S. (2000). Widespread occur- 
rence of structurally diverse 
tetraether membrane lipids: evi- 
dence for the ubiquitous presence 
of low-temperature relatives of 
hyperthermophiles. Proc. Natl. 
Acad. Sci. U.S.A. 97, 14421-14426. 
doi: 10.1073/pnas.97.26.14421 

Schouten, S., Hopmans, E. C, Schefufi, 
E., and Sinninghe Damste, J. S. 
(2002) Distributional variations in 
marine crenarchaeotal membrane 
lipids: a new organic proxy for 
reconstructing ancient sea water 
temperatures? Earth Planet. Sci. Lett. 
204, 265-274. doi: 10.1016/S0012- 
821X(02)00979-2 

Schouten, S., Huguet, C, Hopmans, 
E. C, and Sinninghe Damste, 
J. S. (2007). Improved analyti- 
cal methodology of the TEXgg 
paleothermometry by high 
performance liquid chromatogra- 
phy/ atmospheric pressure chemical 
ionization-mass spectrometry. 
Anal. Chem. 79, 2940-2944. doi: 
10.1021/ac062339v 

Sinninghe Damste, J. S., Hopmans, 
E. C, Schouten, S., van Duin, A. 
C. T., and Geenevasen, J. A. J. 
(2002). Crenarchaeol: the char- 
acteristic glycerol dibiphytanyl 
glycerol tetraether membrane 
lipid of cosmopolitan pelagic 
crenarchaeota. /. Lipid Res. 43, 
1641-1651. doi: 10.1194/jlr. 
M200148-JLR200 

Sinninghe Damste, J. S., Ossebaar, 
J., Abbas, B., Schouten, S., and 
Verschuren, D. (2009). Fluxes and 
distribution of tetraether lipids in 
an equatorial African lake: con- 
straints on the application of the 
TEXgg palaeothermometer and 
BIT index in lacustrine settings. 
Geochim. Cosmochim. Acta 73, 
4232-4249. 

Sun, Q., Chu, G., Liu, M., Xie, M., Li, S., 
Ling, Y, et al. (2011). Distributions 
and temperature dependence of 
branched glycerol dialkyl glycerol 
tetraethers in recent lacustrine sed- 
iments from China and Nepal. 
/. Geophys. Res. 116:G01008. doi: 
10.1029/2010JG001365 

Sun, X., Song, C, Wang, E, and Sun, 
M. (1997). Vegetation history of the 
loess plateau of china during the last 
1000,000 years based on pollen data. 
Q. Int. 37, 25-36. doi: 10.1016/1040- 
6182(96)00008-0 



Tierney, J. E., and Russell, J. M. (2009). 
Distributions of branched GDGTs 
in a tropical lake system: impli- 
cations for lacustrine application 
of the MBT/ CBT paleoproxy. 
Org. Geochem. 40, 1032-1036. 
doi: 10.1016/j.orggeochem.2009. 
04.014 

Tierney, J. E., Schouten, S., Pitcher, 
A., Hopmans, E. C, and Sinninghe 
Damste, J. S. (2011). Core and 
intact polar glycerol dialkyl glycerol 
tetraethers (GDGTs) in Sand 
Pond, Warwick, Rhode Island 
(USA): insights into the origin 
of lacustrine GDGTs. Geochim. 
Cosmochim. Acta 77, 561-581. doi: 
10.1016/j.gca.2011. 10.018 

Walsh, E. M., Ingalls, A. E., and Keil, 
R. G. (2008). Sources and transport 
of terrestrial organic matter in 
Vancouver Island fjords and the 
Vancouver- Washington Margin: 
a multi proxy approach using 
5 13 C (org), lignin phenols, and 
the ether lipid BIT index. Limnol. 
Oceanogr. 53, 1054-1063. doi: 
10.4319/lo.2008.53.3.1054 

Wang, H., Liu, W., Zhang, C. L., 
Wang, Z., Wang, J., Liu, Z., et al. 
(2012). Distribution of glyc- 
erol dialkyl glycerol tetraethers 
(GDGTs) in surface sediments of 
Lake Qinghai and surrounding 
soil. Org. Geochem. 47, 78-87. doi: 
10.1016/j.orggeochem.2012.03.008 

Wang, Y, Cheng, H., Edwards, R. 
L., An, Z., Wu, J., Shen, C.-C, 
et al. (200 1 ) . A high-resolution 
absolute-dated late Pleistocene 
monsoon record from Hulu Cave, 
China. Science 294, 2345-2348. doi: 
10.1126/science.l064618 

Weijers, J. W. H., Bernhardt, B., 
Peterse, R, Werne, J. P., Dungait, 
J. A. J., Schouten, S., et al. 
(2011). Absence of seasonal 
patterns in MBT-CBT indices 
in mid-latitude soils. Geochim. 
Cosmochim. Acta 75, 3179-3190. 
doi: 10.1016/j.gca.2011.03.015 

Weijers, J. W. H., Schouten, S., 
Hopmans, E. C, Geenevasen, J. 
A., David, O. R. P., Coleman, J. 
M., et al. (2006a). Membrane 
lipids of mesophilic anaero- 
bicbacteria thriving in peats have 
typical archaeal traits. Environ. 
Microbiol. 8, 648-657. doi: 
10.1 1 1 l/j.l462-2920.2005.00941.x 

Weijers, J. W. H., Schouten, S., 
Spaargaren, O. C, and Sinninghe 
Damste, J. S. (2006b). Occurrence 
and distribution of tetraether 
membrane lipids in soils: impli- 
cations for the use of the TEXg6 
proxy and the BIT index. Org. 
Geochem. 37, 1680-1693. doi: 
10.1016/j.orggeochem.2006.07.018 



Frontiers in Microbiology | Terrestrial Microbiology 



July 2013 | Volume 4 | Article 199 | 8 



Jia et al. 



Tetraether records in loess-paleosol sequence 



Weijers, J. W. H., Schouten, S., van 
den Donker, J. C, Hopmans, E. 
C, and Sinninghe Damste, J. S. 
(2007). Environmental controls 
on bacterial tetraether membrane 
lipid distribution in soils. Geochim. 
Cosmochim. Acta 71, 703-713. doi: 
10.1016/j.gca.2006.10.003 

Xiao, J., Porter, S. C, An, Z., Kumai, 
H., and Yoshikawa, S. (1995). Grain 
size of quartz as an indicator of 
winter monsoon strength on the 
Loess Plateau of Central China dur- 
ing the last 130,000 yr. Q. Res. 
43, 22-29. doi: 10.1006/qres.l995. 
1003 

Xie, S., Pancost, R. D., Chen, L., 
Evershed, R. P., Yang, H., Zhang, 
K., et al. (2012). Microbial lipid 
records of highly alkaline deposits 
and enhanced aridity associ- 
ated with signifi cant uplift of 
the Tibetan Plateau in the Late 



Miocene. Geology 40, 291-294. doi: 
10.1130/G32570.1 

Zech, R., Gao, L., Tarozo, R., and 
Huang, Y. (2012). Branched glyc- 
erol dialkyl glycerol tetraethers 
in Pleistocene loess-paleosol 
sequences: three case studies. 
Org. Geochem. 53, 38-44. doi: 
10.1016/j.orggeochem.2012.09.005 

Zhang, C. L., Wang, J., Wei, Y, Zhu, 
C, Huang, L., and Dong, H. (2012). 
Production of branched tetraether 
lipids in the lower Pearl River and 
estuary: effects of extraction meth- 
ods and impact on bGDGT prox- 
ies. Front. Microbiol 2:274. doi: 
10.3389/fmicb.2011.00274 

Zhao, Y, Yu, Z., Chen, E, Zhang, J., and 
Yang, B. (2009). Vegetation response 
to Holocene climate change in 
monsoon-influenced region of 
China. Earth Sci. Rev. 97, 242-256. 
doi: 10.1016/j.earscirev.2009.10.007 



Zhou, L. P., Oldfield, E, Wintle, 
A. G., Robinson, S. G., and 
Wang, J. T. (1990). Partly 
pedogenic origin of magnetic 
variations in Chinese loess. Nature 
346, 737-739. doi: 10.1038/ 
346737a0 

Zhu, C, Weijers, J. W. H., Wagner, 
T., Pan, J.-M., Chen, J.-E, and 
Pancost, R. D. (2011). Sources 
and distributions of tetraether 
lipids in surface sediments across a 
large river dominated continental 
margin. Org. Geochem. 42, 376-386. 
doi: 10.1016/j.orggeochem.2011. 
02.002 

Conflict of Interest Statement: The 

authors declare that the research was 
conducted in the absence of any com- 
mercial or financial relationships that 
could be construed as a potential con- 
flict of interest. 



Received: 26 March 2013; paper pend- 
ing published: 29 April 2013; accepted: 
25 June 2013; published online: 15 July 
2013. 

Citation: Jia G, Rao Z, Zhang J, Li Z 
and Chen F (2013) Tetraether biomarker 
records from a loess-paleosol sequence 
in the western Chinese Loess Plateau. 
Front. Microbiol. 4:199. doi: 10.3389/ 
fmicb.2013.00199 

This article was submitted to Frontiers 
in Terrestrial Microbiology, a specialty of 
Frontiers in Microbiology. 
Copyright © 2013 Jia, Rao, Zhang, 
Li and Chen. This is an open-access 
article distributed under the terms of the 
Creative Commons Attribution License, 
which permits use, distribution and 
reproduction in other forums, provided 
the original authors and source are 
credited and subject to any copyright 
notices concerning any third-party 
graphics etc. 



www.frontiersin.org 



July 2013 | Volume 4 | Article 199 | 9 



